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Abstract 



We explore the phase diagram of strongly interacting matter as a func- 
tion of temperature and baryon number density, using a class of models for 
two-flavor QCD in which the interaction between quarks is modelled by that 
induced by instantons. Our treatment allows us to investigate the possible 
simultaneous formation of condensates in the conventional quark-anti-quark 
channel (breaking chiral symmetry) and in a quark-quark channel leading 
to color superconductivity: the spontaneous breaking of color symmetry via 
the formation of quark Cooper pairs. At low temperatures, chiral symmetry 
restoration occurs via a first order transition between a phase with low (or 
zero) baryon density and a high density color superconducting phase. We find 
color superconductivity in the high density phase for temperatures less than 
of order tens to 100 MeV, and find coexisting (qq) and (qq) condensates in 
this phase in the presence of a current quark mass. At high temperatures, 
the chiral phase transition is second order in the chiral limit and is a smooth 
crossover for nonzero current quark mass. A tricritical point separates the 
first order transition at high densities from the second order transition at 
high temperatures. In the presence of a current quark mass this tricritical 
point becomes a second order phase transition with Ising model exponents, 
suggesting that a long correlation length may develop in heavy ion collisions 
in which the phase transition is traversed at the appropriate density. 
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I. INTRODUCTION 



Strongly interacting matter is expected to undergo a phase transition or crossover to a 
quark-gluon plasma phase both at high temperature and at high baryon density. Since QCD 
is asymptotically free, when either the temperature T or the Fermi momentum pp is high the 
effective coupling for typical scattering processes with momentum transfer of order T or pp 
is small, and one therefore expects a new phase of matter in which color is screened rather 
than confined and chiral symmetry is restored. The exploration of the phase diagram is of 
fundamental interest and has applications in cosmology, in the astrophysics of neutron stars 
and in the physics of heavy ion collisions. The zero density, high temperature axis of the 
phase diagram is much better explored than the zero temperature, high density axis, since 
lattice Monte Carlo techniques are well-suited to nonzero temperature but not so well-suited 
to nonzero chemical potential. Recent work suggests a rich phase structure at nonzero 
density. First, chiral symmetry restoration proceeds via a first order phase transition and 
ordinary nuclear matter may be described as being in the mixed phase of this transition 
0. Second, above the transition, the quarks may be weakly coupled but the phase they 
are in is not the trivial one. Even at weak coupling, any attractive quark-quark interaction 
leads to the formation of a condensate of Cooper pairs, which spontaneously breaks color 
symmetry. Color superconductivity has been studied using perturbative methods valid at 
asymptotically high density ||, and in a one-flavor model [f|. The authors of Refs. |],|2| 
have used a model in which the interaction between quarks in two-flavor QCD is modelled 
by that induced by instantons |J to estimate the magnitude of the diquark condensate in 
the phase in which chiral symmetry is restored. We first generalize the variational methods 
of to a formalism in which we derive a bosonic effective action for the degrees of freedom 
which condense.^] For the present investigation, we evaluate the effective action using a mean 
field approximation. Using the method we describe here, we can analyze circumstances in 
which a chiral symmetry breaking (ipip) condensate and a color breaking (ipip) condensate 
coexist. Furthermore, we introduce nonzero temperature and current quark mass. In this 
paper, therefore, we are able to explore the phase diagram as a function of temperature, 
density or chemical potential, and quark mass. 

The phase diagram which we uncover has striking qualitative features, several of which 
we expect to generalize beyond the model which we consider and to provide a good guide 
to the physics of two-flavor QCD. At low temperatures, chiral symmetry restoration occurs 
via a first order phase transition between a phase with low (or zero) baryon density and 
a high density phase with a condensate of quark-quark Cooper pairs in color antitriplet, 
Lorentz scalar, isospin singlet states.0 We find color superconductivity at temperatures 
T < T C A where T C A is of order tens to almost one hundred MeV, and at most a factor of four 
higher if we allow interactions other than that induced by instantons. This suggests that 
color superconductivity may arise in heavy ion collisions, in which the necessary densities 



1 Diquark condensation has also been sought in the vacuum using similar methods ||. 

2 There are also indications [jlj of a color 6, Lorentz axial vector, isospin singlet condensate which 
is many orders of magnitude smaller than the condensates we treat. 
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are likely accompanied by temperatures of order 100 MeV. The transition we find at T C A is 
second order for the present model, but the order of this transition may change once gauge 
field fluctuations are taken into account. We study the competition between chiral and 
superconductor condensates, in particular when the current quark mass is nonzero and both 
coexist. For small quark mass, we find that the color superconductivity is almost unaffected 
by the presence of a chiral condensate. 

In the chiral limit, we find a second order chiral transition at zero chemical potential 
and a first order chiral transition at zero temperature. Consequently, these meet at a 
tricritical point at a particular nonzero temperature T tc and chemical potential /x tc as was 
first discussed in a different phenomenological model in Ref. ]7). We point out that the 
physics in the vicinity of the tricritical point is described by a 6 theory, and the transition 
should therefore be characterized by the mean field exponents which we calculate and by 
logarithmic corrections to scaling ||. If two-flavor QCD has a second order transition 
at high temperatures and a first order transition at high densities, then it will have a 
tricritical point in this universality class. Away from the chiral limit, the second order chiral 
transition turns into a smooth crossover, while the first order transition remains first order. 
Of particular interest is the fact that the tricritical point becomes an Ising second order 
transition. This raises the possibility that even though the pion is massive in nature, long 
correlation lengths may arise in heavy ion collisions which traverse the chiral transition near 
T tc and /i tc . Much of the quantitative physics associated with the first order transition is 
quite model-dependent. However, the physics of the region near the tricritical point is both 
universal and well-described in mean field theory. 

II. MODELS 

Our goal is an exploration of the phase diagram of two flavor QCD in a context which 
allows us to describe likely patterns of symmetry breaking and to make rough quantitative 
estimates. As a tractable model, we consider a class of fermionic models for QCD || at 
nonzero temperature and baryon number density JTO] where the fermions interact via the 
instanton-induced interactions between light quarks || . In the Euclidean functional integral 
for the effective action, nonzero temperature T is implemented via antiperiodic boundary 
conditions for fermionic fields in the 'time' direction with periodicity 1/T. In the presence 
of a nonvanishing chemical potential [i for (net) quark number, the quadratic part of the 
Euclidean action 5* in momentum space readsQ 

with an average current quark mass m. The zeroth component of the momentum q = (q°, q) 
is discrete with q°(n) = (2n+l)7rT for the fermionic Matsubara modes. Indices % — 1, . . . , Nf 
and a — 1, . . . , N c denote flavor and color, respectively. Spinor indices are suppressed and 



3 We use the Euclidean conventions employed in Ref. fll|| . 
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repeated indices are summed. We will write expressions in terms of Nj and N c throughout, 
even though we will only be concerned with Nf = 2 and N c = 3 in this paper.0 

We add a four-fermion interaction to (]T|) which has the color, flavor and Lorentz structure 
of the instanton vertex of two-flavor QCD ||. This interaction properly reflects the chiral 
symmetry of QCD: axial baryon number is broken, while chiral SU(2) l x SU(2)r is respected. 
Color 577(3) is realized as a global symmetry. One could add other four-fermion interactions 
that respect the unbroken symmetries of QCD, including in particular that induced by one 
gluon exchange. We will comment on its effects below. Using the instanton vertex alone 
is only the simplest way of breaking all the symmetries QCD breaks while respecting those 
it respects. We note that with the help of appropriate Fierz transformations the instanton 
interaction can be decomposed into two parts, where one part contains only color singlet 
fermion bilinears and the other part contains only color 3 bilinears: Si = Sj + Sj 3c ^ with 

sj lc) =g 1 tJ2[ S% {-CW<M; -p]o ia )[M;p] - OfofoM; -p\o { ^3\p\ 

= G 2 T J {-Ol^pjO^p} + Ofatf-MOfpyatop]} • (2) 

Here p = (2mrT,p) due to the bosonic nature of the fermion bilinears, and t denotes the 
operation of Euclidean reflection. The bilinears 0( CT ), Of*,, O^i) and O?^ with a = 1,2,3 
carry the quantum numbers associated with the scalar isosinglet (a), the pseudo-scalar 
isotriplet (tt), the pseudo-scalar isosinglet (r/) and the scalar isotriplet (a ), respectively. 
Similarly, the bilinear Ofy (O^L\), with the color index a = 1, 2 or 3, carries the quantum 
numbers of the color antitriplet scalar (pseudo-scalar) diquark. (See (0) below for explicit 
expressions.) The instanton interaction introduces only one coupling. If we take this to 
be Gi, then G2 = G\/(N C — 1) in (|2|). We have generalized the interaction to allow G\ 
and G2 to take on independent values, although we will always assume that they have the 
same sign, as for the instanton interaction. We discuss choices of G\ and G2 at greater 
length below. We note from the signs in Sj that if we choose the sign of G\ such that the 
interaction in the a channel is attractive, so that chiral symmetry breaking is favored, we 
may expect condensation in the 7r and scalar diquark channels also. In the chiral limit, 
one can always make a rotation such that there is no tt condensate. A condensate in the 
pseudoscalar diquark channel would break parity spontaneously, but this seems not to be 
favored by the interaction (|2]). We will use the model to explore condensates in the a and 
scalar diquark channels. 



4 We work throughout with the same chemical potential /i for both up and down quarks. Color 
superconductivity in a system with more down than up quarks, say by about 15% as in a lead 
or gold nucleus, has been considered in Ref. [12]. The superconductor gap is somewhat reduced, 
and this suggests that in a heavy ion collision, the system can lower its energy by equalizing the 
down/up ratio in the densest region of the plasma in order to minimize the gap. This suggests 
that these dense regions will expel negative pions toward the periphery. 12] 
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In order to mimic the effects of asymptotic freedom, the interaction has to decrease 
with increasing momentum. We follow Ref. and implement this via a form factor F(q) 
in momentum space for the three-momenta of each of the fermions, with F(0) = 1 and 
F(q) <C 1 for q 2 ^> A 2 . Here A is some effective QCD cutoff scale which one might anticipate 
to be in the range 1/2 — 1 GeV. With form factors included, the fermion bilinears in the 
above interaction read 

W [^;p] = -iTY, I ^F(q)F(p-q)^ a (-q)^ a (p-q) , 

rf^Z J ( 27r ) 

0? s) [<M=T£ / ^F(-^F(-^+g^(^ T )V-p + g)C T 5 6^%.^ 7 (-?) , (3) 

for the cr and for the scalar diquark and similarly for the other bilinears. C denotes the 
charge conjugation matrix. 



III. THERMODYNAMIC POTENTIAL 

The fermionic effective action F (the generating functional for 1PI Green functions) 
determines the field equations that contain all quantum effects. In thermal and chemical 
equilibrium, T depends on the temperature T and chemical potential /i. Here we are in- 
terested in a computation of the effective action that determines the field equations for the 
expectation values of the fermion bilinears 0^[ij), ip; x] and 0? S M; x] whose Fourier compo- 
nents are given in (H). We denote these expectation values by (ipip) and (tpip), respectively. 
The chiral condensate (ipip) is an order parameter for chiral symmetry breaking for vanishing 
current quark mass m. For nonzero (ipip), the diquark bilinear has a color index a which 
chooses a direction in color space, thus breaking color symmetry SU(3) — ► S77(2).f] In the 
following we compute the effective action as a function of the two (space dependent) order 
parameters.^ Its behavior for constant order parameters yields the effective potential which 
corresponds at its extrema to the thermodynamic potential and encodes the equation of 
state. 

For our purposes, it is advantageous to introduce bosonic collective fields into the func- 



tional integral which defines the fermionic effective action ||13|| . These collective fields carry 
the quantum numbers of the fermion bilinears appearing in (^]). We employ a Hubbard- 
Stratonovich transformation, in which the collective fields are introduced into the functional 
integral by inserting the identities 



5 The £7(1) of electromagnetism is spontaneously broken but there is a linear combination of elec- 
tric charge and color hypercharge under which the condensate is neutral which therefore generates 
an unbroken U(l) gauge symmetry. 

6 We note that the diquark condensate breaks the gauge symmetry and is therefore not a gauge 
invariant order parameter. In this respect the situation resembles that in the electroweak sector. 
As in that case, there is the possibility that a transition associated with a rapid change in physical 
properties may proceed through a smooth crossover without any thermodynamic singularity. 
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(4) 

and similar identities for the other bilinears. (Here N and N are field independent nor- 
malization factors.) The terms quadratic in the fermion bilinears in (f|) cancel the original 
four-fermion interaction @. As a result, the four-fermion interaction is cast into a Yukawa 
interaction between fermions and collective fields and a mass term for the collective fields. 
The fermionic fields then appear only quadratically and can be integrated out exactly, leaving 
an effective action for the bosonic collective fields alone. We apply a saddle point expansion 
to the resulting effective action. In particular, standard mean field results correspond to 
the lowest order in this expansion. Having done the Gaussian integral for the fermions, the 
saddle point effective action reads 

V[<p, A, A*; T,fi] = - tr In H 21 - i tr In [l - (H 21 )- 1 H 22 (H l2 y l H n 



d 3 p 
(2^)3 



1 >(-p)0(p) + ^rK (p)a°(p) 



4Gi r/T ^' 4G 2 



(5) 



where, for constant fields (p(x) = <fi and (taking the color index a = 3) A 3 (x) = Ag(x) = A 
and Aj = A 2 = 0, the H ab are diagonal in momentum space: 



(H n )f(q,q') = -C l5 e a ^F(q) 2 A5 3 (q- q')5 n , n , /2ttT = -(H 22 ^f(q,q') , 
(H 12 )f(q, q') = -[C(-rq, + tf*n)C + i(m + F(g) 2 0)]^^-5 3 (g - q')8 nin , /2*T , 
(if 21 )^(g, ? = [( 7 ^ M + z 7 V) + «(m + F{q) 2 ct>Wn i3 5\q - q')8 n , n , /2-kT . (6) 



The integration over the fermions gives rise to a determinant which has been rewritten as 
the trace of a logarithm in ([5]). The trace involves the momentum integration as well as 
the summation over flavor, color and spinor indices. In principle, the saddle point effective 
action depends on the expectation values of all fermion bilinears occuring in (^). The 
corresponding collective fields would modify the H ab of in a straightforward way and for 
each field a quadratic, mass like, term appears in (|5|). As discussed in the previous section, 
these collective fields are not expected to condense. They therefore only contribute to higher 
order terms in the saddle point expansion. 

In this work, we restrict ourselves to the mean field approximation though the approach 
can be extended to take into account fluctuations around the saddle point. The mean field 
approximation corresponds to the summation of an infinite class of diagrams, but it is not a 
systematic expansion in powers of a small parameter. A nonperturbative approach beyond 
mean field may be performed along the lines of Ref. There, the two flavor effective 



potential has been computed at zero baryon density using renormalization group methods. 
In particular, this approach allows a precise treatment of the effective potential in the vicinity 
of second order fl^ , [TT|| or weak first order [[L5[] transitions, both of which we will encounter 
in Section V. We leave the extension of these methods to nonzero baryon density for future 
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work, as our goal here is to use mean field results as a qualitative guide in exploring the 
phase diagram. 

We evaluate the effective action in the saddle point approximation for constant field 
configurations. This yields the effective potential or thermodynamic potential Q = YT/V. 
A convenient way to extract the effective potential from ([5]) is to consider the derivatives of 
T with respect to <p an d A. After performing the traces and, in particular, the sums over 
Matsubara frequencies, Q can be obtained by integration. We find 



fi(0,A;//,T) 



1 



4Gi 4G 2 
q aq 



A 2 



2N F J |-| (Ar c - 2 ){^ + Tln(l + exp 
o ^ 

+Tln + exp 



(E+-n)/T 
(E^ + n)/T])} 



+ (£ -fJL)y/l + F 4 A 2 /(^ - flf + (E4 + (JL)y/l + F 4 A 2 /(£^ + /i) 2 

+ 2Tln (l + exp 



+ 2Tln 1 + exp 



where F(q) is the form factor and 



E^q) = ^q z + (m + F(q)< 



constant (7) 



(8) 



The constant in (|7|) does not depend on /i and T and is chosen such that the pressure of the 
physical vacuum is zero. Extremizing Q leads to coupled gap equations for the two order 
parameters, whose solutions we call (f> and A . That is, 



;A=A n 



dA 



= 



(9) 



A ;A=A 



where O an d A are related to the chiral condensate and the condensate of Cooper pairs by 

0o = 2Gi(^), A = 2G 2 (^)- (10) 



At its extrema, the potential is related to the energy density e, the entropy density s, the 
quark number density n and the pressure P by 



fi(0O) Aq; fJ>,T) = e — Ts — fin = —P . 



(11) 



Unless otherwise stated, the results we quote in our exploration of the phase diagram are 
obtained using the smooth form factor 



F(q)=A 2 /(f + A 2 ) 



(12) 



with A = 0.8 GeV and with G a A 2 = 6.47 fixed by requiring O = 0jf c = 0.4 GeV at 
fj, = T = in order to obtain a reasonable, albeit qualitative, phenomenology. For most of 



6 







FIG. 1. The thermodynamic potential O (in GeV 4 ) as a function of <p an d A at T = and 
/i = 0.292 GeV. The two degenerate minima have 4> = <^f c = 0.4 GeV, A = and <p = 0, A = 0.072 
GeV. 

the discussion we will consider G2 = 3Gi/4, which will be motivated at length below. We 
note that the qualitative features which we address do not depend on this specific choice 
and we discuss the effects of different choices below. 

Although we describe our results in the next sections, it is worth pausing here to em- 
phasize how we use Q. As T or fi or m changes, Q can have several local minima in the 
(<f), A) plane. The lowest minimum describes the lowest free energy state and is favored. As 
an example, which we discuss at length in the next section, in Figure 1 we plot fl(4>, A) at 
T = and fi = 0.292 GeV. One observes two degenerate minima corresponding to a first 
order phase transition at which two phases have equal pressure and can coexist. We discuss 
the symmetry breaking patterns, their effects and the varied phase transitions which occur 
as a function of temperature and chemical potential or density in Sections IV- VI. 

Before proceeding, we return to our postponed discussion of the appropriate choice for 
G2/G1, having fixed G\ phenomenologically. If we take G2 to be independent of G\, there 
is a phenomenological upper bound which G2 must satisfy. It is only for G2 < 2G\ that 
the minimum of the potential at \i — T — has A = 0. For larger G2, color would be 
spontaneously broken in the vacuum (see also [§]). (As G2 increases above 2G\, the vacuum 
state moves continuously in the (</>,A) plane from the A = axis to the = axis.) It 
is not unreasonable to leave G2 as a free parameter satisfying G2 < 2Gi for the following 
reason. One could extend our model by introducing four fermion interactions different from 
(0) with new coupling constants, as long as the correct symmetries are respected. These 
could again be Fierz transformed into products of fermion bilinears and treated as described 
above. Assuming that condensation occurs only in the a and scalar diquark channels, we 
would end up with a mean field thermodynamic potential precisely of the form (|7]), but with 
G\ and G2 replaced by rescaled couplings G[ and G' 2 . However, there is a caveat if a vector 
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interaction of the form (ip^-ip) 2 is present. As with any four fermion interaction, this can 
rescale Gi and G 2 . However, because (ipj ip) = n ^ 0, there is an additional qualitative 
effect. In mean field, the vector interaction translates into an n 2 term added to ([!]) and an 
effective shift in the chemical potential proportional to the density n |I(J. We will address 
the consequences of such an interaction in Section IV. 

Suppose that, instead of leaving GifG\ unspecified, we wish to treat the instanton inter- 
action for which G^ — Gxj (N c — 1). If we could analyze the model without approximation, 
we would use this ratio. However, the decomposition of the instanton interaction in (0) is, 
of course, not unique. It can easily be Fierz transformed at the cost of generating products 
of other bilinears which are, for example, in color octet or 6 representations. Though this 
is of no concern were we to treat the model exactly, it introduces an ambiguity in the mean 
field treatment we employ. This is because in addition to generating new bilinears, the Fierz 
transformation changes the couplings in front of the bilinears in (fj). In complete analogy to 
the above discussion of new interactions, this would change the result (|7|) for Q by a rescaling 
of Gi and G 2 . This ambiguity is a familiar one in studies of chiral symmetry breaking in 
NJL models, where it is well known that a mean field calculation of the chiral condensate 
reproduces the results of summing direct and exchange diagrams if one adds to the original 
interaction a suitably Fierz transformed interaction [M. Applied to (0), this scales G<i to 
zero, removing the diquark bilinears, and introduces color octet terms. This procedure is 
correct only if no diquark condensate forms, as is of course the case at low densities. Because 
we are considering the possibility of diquark condensation in addition to chiral condensation, 
we must then Fierz transform the interaction (including the color octet terms) yet again, this 
time into a form with G\ scaled to zero, containing no quark-ant i- quark bilinears. Adding 
this to the interaction without diquark bilinears,[] we obtain ([|) with 

N 

G 2 = G x — , (13) 

and with color octet and 6 terms which we now discard because they are assumed not to 
condense. This procedure yields gap equations which agree, when either or A is set to zero, 
with those obtained previously pj via extremizing a variational wave function and we use 



the ratio fll3D in all the results we quote unless stated otherwise. The variational method has 
the advantage of avoiding these ambiguities in the mean field approach we are using here. 
Our present analysis has the virtue of being tractable at nonzero temperature, for nonzero 
quark mass and, perhaps most important, in the presence of simultaneous condensation in 
the a and scalar diquark channels. 

IV. ZERO TEMPERATURE PHYSICS 

We begin our analysis of the thermodynamic potential (|7|) at zero temperature and in 
the chiral limit. fl(<p, A; jA,T) is a function of two order parameters. As \i changes, Q can 



7 This is equivalent to the observation that whereas chiral condensation uses the instanton vertex 
to turn a ulur into a (IrcIl, condensation in the scalar diquark channel uses it to turn ulcIl into 
urcIr, and so one must add the exchange interaction in this sense also. 
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FIG. 2. The zero temperature thermodynamic potential f2 (in GeV 4 ) as a function of at 
A = (left panel) and as a function of A at = (right panel) for several chemical potentials. 
The curves correspond to (top to bottom) fi = 0, 0.292, 0.35 GeV. The curves at fi = fio = 0.292 
GeV are sections of Figure 1. 

have several local minima in the (0, A) plane. An example has been given in Figure 1, where 
for fi = fiQ = 0.292 GeV the potential has two degenerate minima. To see how the potential 
changes with fi, in Figure 2 we display sections of the (0, A) plane for A = and = 
respectively. Let us first focus on the left panel. At fi — 0, Q has a minimum at = 0.4 
GeV. As fi is increased, Q at this minimum is unchanged. Since 

n = — — — at fixed T, = O , A = A , (14) 

Ofi 

we conclude that the phase described by this minimum has baryon density zero, as we expect 
for the vacuum. As fi is increased sufficiently, a second local minimum of Q develops with 
= and A / 0. The system remains in the vacuum state with density zero and 0o 
unchanged until the chemical potential is increased to fi = fio where both minima become 
degenerate with Q = —P = as seen in Figure 1. There is now a zero pressure chirally 
symmetric color superconductor phase whose density we denote no which has the same 
fi and P as the vacuum phase. Before fi can be increased any further, the system must 
undergo a first order phase transition during which the density increases from zero to n . 
At all intermediate densities, the system is in a mixed phase in which there are regions with 
density Uq and regions of vacuum. Once the transition is complete, \x can increase further, 
and the lowest Q configuration becomes one with Q < 0, that is with positive pressure, in 
which n > uq. In the top panels of Figure 3, we show the behavior of the global minimum 
(0o, Ao) of Q and the density as a function of /i. We have verified that there are no other 
minima, for example with both <p and A nonzero, as long as we stay in the chiral limit. The 
locations of the extrema of the curves in Figure 2 agree with those obtained in Ref. 0] and 



described also in Ref. [he]. Treating the full problem, one observes that the first order phase 
transition does in fact occur between the two minima in Figure 1, at fi — /io = 0.292 GeV. 

There is an interesting interpretation of this zero temperature first order phase transition 
|I1 JT6|JT7[| . For reasonable choices of parameters, n is greater than nuclear matter density 
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FIG. 3. Left panels show 0o and Ao (namely and A at the global minimum of £1) as functions 
of \i for T = (above) and T = 0.03 GeV (below). The right panels show the quark number density 
n for the same parameters. We describe the results at nonzero temperature in the next section. 



and comparable to the baryon density in a nucleon^ Ordinary nuclear matter is then in the 
mixed phase of this transition and consists of nucleon sized droplets within which n = and 
= 0, surrounded by regions of vacuum. Although it is nice to see that the model forces 
quarks to be located within droplets within which they have zero mass, the model must 
obviously be extended in several respects. First, a short range repulsion between droplets 



must be included, so that small droplets are favored over bigger ones |Tc|. Second, color 
must be gauged if the model is to yield droplets which are color singlets. It is crucial to 
this interpretation that /xo < 0o ac - m the chiral limit, 0o corresponds to the constituent 
quark mass and therefore if \i is increased above 0™ the (formerly) vacuum phase becomes 
populated with a gas of constituent quarks. If /i < 0o ac , the first order transition to the high 
density phase which occurs at /i = fi preempts this. Within the present model, /i < 0o ac 



If we identify no with the quark number density in a nucleon and take this to be three times 

1/3 

that in nuclear matter, we should require n = 0.22 GeV. The choice of parameters we are using 

1/3 

for illustration yields a smaller value n Q = 0.18 GeV in the chiral limit, as seen in Figure 3. We 

1/3 

shall see in Section VI that we obtain n = 0.22 GeV for a quark mass of 10 MeV. 
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only for A < 2.2 GeV. Furthermore, adding a four fermion interaction to model one gluon 
exchange introduces an effective shift in \i as discussed in Section III, and this increases 



the value of /io pi. With vector interactions included, one can construct models with 



reasonable values of a suitably defined A in which /i > 0™° 113 • m sucn a model, the 



first order transition occurs between a gas of constituent quarks with a nonzero density n_ 
and a gas of massless quarks with a density n + > n_, and n_ can be greater than nuclear 
matter density for reasonable choices of parameters. If one then gauges color, the gas of 
constituent quarks which the model describes at densities below n_ is expected to become 
a gas of baryons. We now turn to nonzero temperature and one of the things which we will 
see is that models with n_ = and with n_ > are less different than they appear to be at 
zero temperature. 



V. THE PHASE DIAGRAM 

In Figure 4, we present the phase diagram of the model as a function of \i and T, and 
as a function of n and T. We devote this section to a discussion of the physics shown in 
these figures. Let us first focus on Figure 4a, and ignore the quark number density n. We 
find that at /i = 0, the chiral phase transition as a function of increasing temperature is 
second order. In our present mean field treatment, it is of course a mean field transition. 
A second order chiral phase transition in two flavor QCD at nonzero temperature has the 
appropriate symmetry to be in the 0(4) universality class p!8|-[2T|], and a treatment of the 



present model which went beyond mean field theory, like that in Ref . JTTJ] , would find critical 
exponents characteristic of this universality class. There is no reason to expect a small 
chemical potential to change this result, since this introduces no new massless degrees of 
freedom in the effective three dimensional theory which describes the long wavelength modes 
near T = T c . This argument has been made at greater length in Ref. [p2| |. However, it is 



possible that as a parameter in the theory is changed, the quartic coupling in the effective 
three dimensional theory can become negative, making the transition first order. Where 
the quartic coupling is zero, one has a tri critical point. (This is precisely what is expected 
to occur at /i = at a particular value of the strange quark mass. |I9|- |21|J23[] ) Since we 



have seen that the zero temperature, nonzero \i transition is first order and we have found a 
second order transition at high temperature, we expect and find a tricritical point where the 
second order transition becomes first order. The tricritical point occurs at fi = fi tc = 211 
MeV and T = T tc = 101 MeV for the parameters we are using. 

We have verified that at the tricritical point Q ~ <p 6 as we expect, based on the arguments 
above.0 In the vicinity of the tricritical point, Q therefore has the form 

fi(0, 0; /i, T) = fi(0, 0; A i t0 , T tc ) + + + - H , (15) 

where the coefficient h of the linear term is proportional to the current quark mass. The 
coefficients a and b are both zero at the tricritical point, and barring accidental cancellations 



9 In fact, tt(<f>, A = 0; fj, tc , T tc ) = O.O330 6 - 0.000392, all in GeV. 
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FIG. 4. Phase diagram as a function of T and /i, and as a function of T and n 1 / 3 (all in GeV). 
This section can be viewed as one long caption for these figures. The solid curves are second order 
phase transitions; the dashed curves describe the first order transition. 

both will be linear in both (T — T tc ) and (p — /x tc ). The /i and T dependence of c is not 
important, as c does not vanish at the tricritical point, and it is convenient to set c = 1. It 
is easy to verify that, near the tricritical point, the line of second order transitions is given 
by a = 0, b > and the line of first order transitions is given by a = 36 2 /16, b < 0. Minima 
of Q are described by the scaling form 



h 




(16) 



From this, we read off the exponents 6 — 5, 1//3 — 4, and 4>t//3 = 2 where <p t is called 
the crossover exponent because tricritical (as opposed to first or second order) scaling is 
observed for b < or*. For more details see Refs. |]8|j2T[|. The relation between a and b and 
(T — T tc ) and (fj, — /it c ) depends on the slope of the line of phase transitions which end at 
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the tricritical point in Figure 4a and is not universal. In general, our model should not 
be used quantitatively. However, if we trust the qualitative feature that the transition is 
second order at high temperatures and first order at low temperatures, then QCD with two 
flavors will have a tricritical point in the same universality class as that in our model. Since 
d — 3 is the critical dimension for <p G theory, we expect that a complete treatment would 
yield mean field exponents in agreement to those in our model, although there would be 
logarithmic corrections to the scaling law ( |TE| ) || . The reason we have described the physics 
of the tricritical point carefully is that, because we have a 6 theory in three dimensions, 
the universal critical exponents will be given quantitatively by those we find in mean field 
theory. 

At low temperatures, one has color superconductivity at chemical potentials immediately 
above the first order phase transition. However, the superconducting order parameter A is 
reduced by nonzero temperature, as shown in Figure 3. In Figure 4 we see that for the 
parameters we are using, color superconductivity persists up to temperatures as high as 45 
MeV for the most favorable chemical potential. The energy gap in the color superconductor 
phase (that is, the energy to create a quasiparticle-quasihole pair) is 2A F(fi) 2 . BCS theory 
predicts that, at weak coupling, the critical temperature T A at which the superconductivity 
vanishes should be given by 0.57F(/i) 2 A (T = 0). Our results are within a few percent of 
this. This suggests that T A is sensitive to the shape of the form factor. For example, we 
have verified that if we use a form factor with a sharp cutoff we can easily increase T A by a 
factor of two, although this tends to increase T c for the finite temperature transition at zero 
chemical potential also. If we generalize beyond the instanton interaction, and therefore do 
not require but simply impose G2 < 2G± as required if color SU (3) is to be a symmetry 
of the vacuum, then we find that values of A and T C A can arise which are up to four times 
larger than those we show in our figures. 

The transition at nonzero temperature at which color superconductivity is lost is second 
order in our model. Because of the shape of the T A curve as a function of fi, there are 
temperatures (like T = 30 MeV in Figure 3) for which, as fi is increased, one first completes 
the first order chiral transition, then has a region with <p = A = 0, then has a second 
order transition to a color superconducting phase, before finally having another second order 
transition to a phase with no condensates. If one neglects gluons, one could study this second 
order transition by doing a renormalization group analysis to first determine whether there 
is a second order transition as in our mean field analysis, and then to determine the universal 
exponents due to the long wavelength fluctuations of A, associated with the breaking of the 
global color symmetry. We expect, however, that the second order nature of the transition 
at T = T A can change once gauge fields are introduced. For example, the electroweak phase 
transition is second order when only the Higgs field is included, but can be either first order 
or a smooth crossover once gauge fields are taken into consideration. We leave a Ginzburg- 
Landau treatment of the physics near T A (including both gluons and the superconductor 
order parameter A) to future work. 

By not showing n, Figure 4a obscures important features seen in Figure 4b, including 
for example the fact (as discussed above) that at zero temperature, one of the phases in 
equilibrium at the first order transition has zero density. At any nonzero temperature 
which is low enough that the transition is first order, the transition occurs between a gas of 
constituent quarks with density n_ and a gas of massless quarks with density n + . In other 
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FIG. 5. As for Figure 3 at zero temperature, but with a current quark mass m = 10 MeV. We 
observe coexisting cpQ ~ (ipip) and Aq ~ {ipip) condensates in the high density phase. 



words, in the mixed phase one has droplets containing massless quarks as before, but they 
are now surrounded by a dilute gas of constituent quarks, rather than by vacuum. Once 
gluons are introduced, both the droplets of the n = n + phase and the constituent quarks in 
the low density phase are expected to form color singlets. Thus, the description of baryons 
as droplets is less complete at nonzero temperature than at zero temperature. 

It is clear that there is no model-independent significance to the dashed lines in Figure 
4b which mark the boundaries of the mixed phase. In nature, with gluons and with repul- 
sive short range interactions between nucleons, the only curve which can be defined is the 
one which bounds the region in which chiral symmetry is broken. Locating this boundary 
requires further analysis of the first order transition, but it seems likely jljl6,24]| that it is 
associated with the percolation of the droplets of the high density phase. More precisely, 
once one no longer has regions of the low density phase (with nonzero 0o) which are of infi- 
nite extent, chiral symmetry will be restored. Thus, the density at which chiral symmetry 
is restored must be somewhat to the left of n + . Similarly, color superconductivity sets in 
at densities above which there are regions of the high density phase of infinite spatial ex- 
tent. Chiral symmetry restoration is associated with the "unper eolation" of the low density 
phase, while the onset of color superconductivity is associated with the percolation of the 
dense phase, and both occur at densities somewhat below n + . None of this physics is visible 
in Figure 4a, since it all occurs within the mixed phase region. It is important to realize, 
therefore, that if one is interested in the experimentally relevant question of whether or not 
a first order transition occurs as n is increased above nuclear matter density, this analysis 
does not provide an answer. The percolation transition could be smooth or second order, 
even though it is occuring "within" the first order phase transition described by the model. 
What is robust, and as we shall see in the next section may be quite significant, is that 
regardless of how the physics as a function of n works within the mixed phase region, there 
is a tricritical point in the phase diagram. 
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VI. PHYSICS AWAY FROM THE CHIRAL LIMIT 



In the previous section, we have explored the phase diagram at nonzero temperature 
and baryon density in the chiral limit. We now turn on a small but nonzero quark mass 
m, as in nature. We see from Figure 5 that in the low temperature phase, the quark mass 
affects 0o, but leaves A = 0. We observe that 0Q ac , /xo and the critical temperature for 
the chiral transition at zero density all increase with m. One might reasonably choose to 
reduce G\ as m is increased in order to keep 0™° fixed. With G 2 jG\ fixed, this would 
reduce G 2 and would therefore reduce A . We do not do this because our goal is to find 
any effects of the presence of the chiral condensate on the superconductor condensate and 
we do not wish these effects to be entangled with the effects of a reduction in G 2 . In 
the phase with density no = (0.217 GeV) 3 attained immediately upon completion of the 
phase transition, both condensates have comparable magnitudes with O = 0.136 GeV and 
A = 0.110 GeV. We observe only a very small dependence of A on m, as long as the 
constituent quark mass m + O is smaller than u. (Compare Figure 5 to Figure 3, where 
m = 0.) This insensitivity can be understood by noting that although the Cooper pairs 
have low momentum, they are formed from quarks which have momenta close to the Fermi 
surface. Adding a quark mass m u does not significantly affect the density of states or 
the interactions of the quasiparticles with momenta of order «. This is consistent with our 
finding that the temperature T A at which the superconductivity is lost is almost unchanged. 

The phase transition at T = T A remains second order in the presence of a nonzero quark 
mass. This result can be understood by expanding the effective potential at small and A. 
Whereas a quark mass adds a term to the effective potential which is linear in at small 
0, it does not introduce a term linear in A. It can modify the coefficient of the A 2 term, 
which changes the critical temperature for the transition but leaves it second order. 

To conclude our discussion of color superconductivity, we remark that the phenomenon 
seems quite robust. It is not significantly disturbed by the simultaneous presence of a 
chiral condensate as long as the constituent quark mass is less than the chemical potential. 
Although the phenomenon is robust, our estimate of the critical temperature at which it 
vanishes is less so. For the parameters we have used for illustrative purposes, T A is 30-45 
MeV, but changes in the form factor can double this and phenomenologically acceptable 
changes in the coupling constant G 2 can multiply T C A by up to a factor of four. A definitive 
answer of whether or not color superconductivity can arise in heavy ion collisions, in which 
the necessary densities are likely accompanied by temperatures of order 100 MeV, must 
await a less qualitative treatment. 

Our model reveals the possibility of a tricritical point in the phase diagram for two-flavor 
QCD in the chiral limit. Physics in the vicinity of the tricritical point is described by a 6 
field theory, whose critical exponents are given quantitatively by the mean field analysis in 
this paper. What, then, happens in this region of temperature and chemical potential in 
the presence of a small quark mass? It is well known that a nonzero quark mass does have 
a qualitative effect on the second order transition which occurs at temperatures above T tc . 
This 0(4) transition becomes a smooth crossover, because terms linear in now arise in Q. A 
small quark mass cannot eliminate the first order transition below T tc . Therefore, whereas 
we previously had a line of first order transitions and a line of second order transitions 
meeting at a tricritical point, with m^Owe now have a line of first order transitions ending 
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at an ordinary critical point, as in the liquid-gas phase diagram. At this critical point, one 
degree of freedom (that associated with the magnitude of <fi) becomes massless, while the 
pion degrees of freedom are massive since chiral symmetry is explicitly broken. Therefore, 
this transition is in the same universality class as the three-dimensional Ising model. 

From many studies of QCD at nonzero temperature, we are familiar with the possibility 
of a second order transition, with infinite correlation lengths, in an unphysical world in which 
there are two massless quarks. It is exciting to realize that if the finite density transition 
is first order at zero temperature, as in the models we have considered, then there is a 
tricritical point in the chiral limit which becomes an Ising second order phase transition 
in a world with chiral symmetry explicitly broken. In a sufficiently energetic heavy ion 
collision, one may create conditions in approximate local thermal equilibrium in the phase 
in which spontaneous chiral symmetry breaking is lost. Depending on the initial density and 
temperature, when this plasma expands and cools it will traverse the phase transition at 
different points in the (//, T) plane. Our results suggest that in heavy ion collisions in which 
the chiral symmetry breaking transition is traversed at baryon densities which are not too 
high and not too low, a very long correlation length in the a channel and critical slowing 
down may be manifest even though the pion is massive. 
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